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Abstract 

We calculate the relic density of the lightest neutralino in a supersymmetric seesaw type-II 
("triplet seesaw") model with minimal supergravity boundary conditions at the GUT scale. The 
presence of a triplet below the GUT scale, required to explain measured neutrino data in this 
setup, leads to a characteristic deformation of the sparticle spectrum with respect to the pure 
mSugra expectations, affecting the calculated relic dark matter (DM) density. We discuss how the 
DM allowed regions in the (mo,Mi/ 2 ) plane change as a function of the (type-II) seesaw scale. 
We also compare the constraints imposed on the models parameter space form upper limits on 
lepton flavour violating (LFV) decays to those imposed by DM. Finally, we briefly comment on 
uncertainties in the calculation of the relic neutralino density due to uncertainties in the measured 
top and bottom masses. 
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I. INTRODUCTION 



Standard cosmology requires the existence of a non-baryonic dark matter (DM) contri- 
bution to the total energy budget of the universe [H, E|- In the past few years estimates of 
the DM abundance have become increasingly precise. Indeed, the Particle Data Group now 
quotes at 1 a c.l. (3| 

tt DM h 2 = 0.105 ±0.008. (1) 

Since the data from the WMAP satellite (3, and large scale structure formation 0| is 
best fitted if the DM is cold, weakly interacting mass particles (WIMP) are currently the 
preferred explanation. While there is certainly no shortage of WIMP candidates (lists can be 
found in many reviews, see for example HQ Eli), the literature i. completely dominated 
by studies of the lightest neutralino. 

Neutrino oscillation experiments have shown that neutrinos have non-zero mass and mix- 



les 0, y, 
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13j and the most recent global fits to all data [lj] confirm again 



ing ang 

that the mixing angles are surprisingly close to the so-called tri-bimaximal mixing (TBM) 
values 



151 ]. In the minimal supersymmetric extension of the standard model (MSSM) with 



conserved R-parity neutrino masses are zero for the same reasons as in the SM. However, it 
was shown long ago that if neutrinos are Majorana particles, their mass is described by a 
unique dimension-5 operator 
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= ~{HL){HL). 



(2) 



All (Majorana) neutrino mass models reduce to this operator at low energies. If / is a 
coefficient 0(1), current neutrino data indicates A < O(10 15 ) GeV. This is the essence of the 
"seesaw" mechanism. There are three different tree-level realizations of the seesaw, classified 
as type-I, type-II and type-IH in [13]. Type-I is the well-known case of the exchange of a 
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leavy fermionic singlet [13,113,120]. Type-II corresponds to the exchange of a scalar triplet 



221] ■ One could also add one (or more) fermionic triplets to the field content of the SM 
23l ] . This is called seesaw type-Ill in ET 



Neutrino experiments at low energies measure only / Ql g/A, thus observables outside the 
neutrino sector will ultimately be needed to learn about the origin of eq. <^j. Augmenting 
the SM with a high-scale seesaw mechanism does not lead to any conceivable phenomenol- 
ogy apart from neutrino masses, but if weak scale supersymmetry exists indirect probes into 
the high energy world might be possible. Two kind of measurements containing such indi- 
rect information exist in principle, lepton flavour violating (LFV) observables and sparticle 
masses. 

Assuming complete flavour blindness in the soft supersymmetry breaking parameters at 
some large scale, the neutrino Yukawa matrices will, in general, lead to non-zero flavour 
violating entries in the slepton mass matrices, if the seesaw scale is lower than the scale at 
which SUSY is broken. This was first pointed out in [24j. The resulting LFV processes have 
been studied in many publications, for low-energy observables such as \i — > ej and /i — e 



conversion in seesaw type-I see for example [25J, 120, 123, 1221, 120, 120, ]31fl, for seesaw type-II 
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32l . l33l |. LFV collider observables have also been studied in a number of papers, see for 



example [33l . 
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Mass measurements in the sparticle sector will not only be necessary to learn about the 
mechanism of SUSY breaking in general, but might also reveal indications about the scale 
of the seesaw mechanism. However, very precise knowled ge o f masses will be necessary 
before one can learn about the high scale parameters [46|, l43|. Especially interesting in 
this context is the observation that from the different soft scalar and gaugino masses one 
can define certain combinations ("invariants") which are nearly constant over large parts of 
mSugra space. Adding a seesaw mechanism of type-II or type-Ill these invariants change in 
a characteristic way as a function of the seesaw scale and are thus especially suited to extract 
information about the high energy parameters |48j|. Note, however, that the "invariants" are 
constants in mSugra space only in leading order and that quantitatively important 2-loop 
corrections exist 
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In this paper we study neutralino dark matter [49l . [50j, iMfl within a supersymmetric 
type-II seesaw model with mSugra boundary conditions. For definiteness, the model we 
consider consists of the MSSM particle spectrum to which we add a single pair of 15- and 
15-plets. This is the simplest supersymmetric type-II setup, which allows one to maintain 
gauge coupling unification [32| and explain measured neutrino oscillation data. 

In mSugra - assuming a standard thermal history of the early universe 1 - only four 
very specific regions in parameter space can correctly explain the most recent WMAP data 
These are (i) the bulk region; (ii) the co-annihilation line; (iii) the "focus point" line 
and (iv) the "higgs funnel" region. In the bulk region there are no specific relations among 
the sparticle masses. However, all sparticles are rather light in this region, so it is already 



53 . In the co-annihilation line 



very constrained from the view point of low-energy data 
the lightest scalar tau is nearly degenerate with the lightest neutralino, thus reducing the 

54j. In the "focus point" line 



neutralino relic density with respect to naive expectations [50 

54, El]] & x oh 2 is small enough to explain Q^Mh 2 due to a rather small value of /i leading to 



an enhanced higgsino component in the lightest neutralino and thus an enhanced coupling to 
the Z° boson. Lastly, at large tan j3 an s-channel resonance pair annihilation of neutralinos 
through the CP-odd higgs boson can become important. This is called the "higgs funnel" 
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region 

The addition of the 15 and 15 pair at the high scale does not, in general, lead to the 
appearance of new allowed regions. However, the deformed sparticle spectrum with respect 
to mSugra expectations leads to characteristic changes in the allowed regions as a function 
of the unknown seesaw scale. We discuss these changes in detail and compare the results to 
other indirect constraints, namely, the observed neutrino masses and upper limits on LFV 
processes. We concentrate on the seesaw type-II scheme, since for mSugra + seesaw type-I 
the changes in the DM allowed regions with respect to pure mSugra are, in general, expected 
to be tiny. 2 



1 In models with non-standard thermal history the relation between sparticle masses and relic density can 
be lost completely 

2 We have confirmed this general expectation with some sample calculations. However, an exceptional case 
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The rest of this paper is organized as follows. In the next section we briefly summarize 
the main ingredients of the model and give a short discussion of mSugra and the expected 
changes in sparticle masses in our setup with respect to mSugra. In section (IHIf l we present 
our numerical results. This is the main section of the current paper, where we discuss in 
detail how the introduction of a 15 changes the predicted DM abundance as a function 
of the seesaw scale. We also confront the DM allowed regions with constraints from non- 
observation of LFV processes and briefly comment on DM in mSugra with a seesaw type-I. 
We then close with a short summarizing discussion in section (TTVl) . 



II. SETUP: MSUGRA AND SU(5) MOTIVATED TYPE-II SEESAW 



In this section we summarize the main features of the model we will use in the numerical 
calculation. We will always refer to minimal Supergravity (mSugra) as the "standard" against 
which we compare all our results. The model consists in extending the MSSM particle 
spectrum by a pair of 15 and 15. It is the minimal supersymmetric seesaw type-II model 
which maintains gauge coupling unification 
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mSugra is specified by 4 continuous and one discrete parameter [57]]. These are usually 
chosen to be m , the common scalar mass, Ml/ 2 , the gaugino mass parameter, A , the 
common trilinear parameter, tan/3 = ^ and the sign of fi. m , My 2 and A are defined at 
the GUT scale, the RGEs are known at the 2-loop level 58 . 

Under SU(3) x SU L (2) x U(1) Y the 15 decomposes as 

15 = S + T + Z (3) 

S ~ (6,1,-^), T~ (1,3,1), Z~(3,2,K 
3 o 

The 577(5) invariant superpotential reads as 

W = -^Y 15 5 • 15 ■ 5 + -=Ai5 H ■ 15 • l H + ^A 2 5^ • 15 ■ 5 H + Y 5 10 • 5 • l H (4) 
v2 \/2 \/2 

+ YiqIO ■ 10 • 5 H + M 15 15 • 15 + M 5 5 

Here, 5 = (d c ,L), 10 = (w c ,e c ,Q), 5 H = (t,H 2 ) and 5 H = (£,#1). Below the GUT scale in 
the SU (5)-broken phase the potential contains the terms 

{Y T LT X L + Y s d c Sd c ) + Y z d c ZL + Y d d c QH x + Y u u c QH 2 + Y e e c LH 1 (5) 

V2 

+ ^(A^Titfx + \ 2 H 2 T 2 H 2 ) + M T T X T 2 + M Z Z X Z 2 + M S S 1 S 2 + fiH^ 
v2 

Yd, Y u and Y e generate quark and charged lepton masses in the usual manner. In addition 
there are the matrices Yt, Ys and Yz. For the case of a complete 15, apart from calculable 
threshold corrections, Y T = Y S = Y z and M T , Ms and M z are determined from M 15 by the 



has been presented recently in [56j, see the more detailed discussion in section ffllljl . 
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RGEs. As long as Mz ~ Ms ~ Mt ~ M i5 gauge coupling unification will be maintained. 
The equality need not be exact for successful unification. 

The triplet T\ has the correct quantum numbers to generate neutrino masses via the 
first term in eq. ([5|). Integrating out the heavy triplets at their mass scale a dimension-5 
operator of the form eq. (J2j) is generated and after electro-weak symmetry breaking the 
resulting neutrino mass matrix can be written as 



2 M T 



(6) 



Here v 2 is the vacuum expectation value of Higgs doublet H 2 and we use the convention 
(Hi) = m u can be diagonalized in the standard way with a unitary matrix U, containing 

in general 3 angles and 3 phases. Note that Y T = U T ■ Y T ■ U is diagonalized by the 
same matrix as m v . This means that if all neutrino eigenvalues, angles and phases were 
known, would be completely fixed up to an overall constant, which can be written as 
10 15 GeV ( 05 cV ) . Thus, current neutrino data requires M T to be lower than the 



A 2 



GUT scale by (at least) an order or magnitude. 

The full set of RGEs for the 15 + 15 can be found in 



32 



and in the numerical calcu- 



lation, presented in the next section, we solve the exact RGEs. However, for a qualitative 
understanding of the results, the following approximative solutions are quite helpful. 
For the gaugino masses one finds in leading order 



MiimsusY) = 7tj v Mi/2- 

a{M G ) 



(7) 



Eq. fZ)) implies that the ratio M 2 /Mi, which is measured at low-energies, has the usual 
mSugra value, but the relationship to Myz is changed. Neglecting the Yukawa couplings 
Y15 (see below), for the soft mass parameters of the first two generations one gets 



n 



^ {MtU f i + fi\M? /2 , 



t=i 




1 + 



a(M G 



An 



h log 



mi 



h + Abi 



1 - 



' a(M G ) n . . Ml 
1+ . (bi + A6«)log 



47T 



(8) 



(9) 



5 



The various coefficients q can be found in [33JJ. The gauge couplings are given as 



3 cos 2 W sin 9 W 



®i{rn SUSY ) 
oti{M T ) 
<Xi{M G ) 



oii{m z ) 



4-rr 2 & 



A/j, 



■■"SUSY 

at (M T ) 

1 - 2^(6, + Afc) log $ 



with fef M and bf 1 being the usual standard model and MSSM coefhcients. A6j = 7 for 
all z in case of a complete 15-plet. 

We can estimate the soft mass parameters given the above formulas for a given choice of 
mo, Mi/2 and M 15 = Mt- We show some arbitrarily chosen examples in fig. (pQ). Note that 
the result shown is approximate, since we are (a) using the leading log approximation and (b) 
two loop effects are numerically important, especially for mg, but not included. The figure 
serves to show that for any M 15 < Mqut the resulting mass parameters are always smaller 
than the mSugra expectations for the same choice of initial parameters (m ,Mi/ 2 ). While 
the exact values depend on (m , M1/2) and on the other mSugra parameters, this feature is 
quite generally true in all of the (m ,Mi/ 2 ) plane. Note, that the running is different for 
the different scalar mass parameters, but the ratio of the gaugino mass parameters Mi/M 2 
always stays close to the mSugra expectation, Mi ~ ^t&n 2 9 wM 2 . 




10 13 10 14 10 15 10 16 10 13 10 14 10 15 10 16 

Mis = M T [GeV] M 15 = M T [GeV] 



FIG. 1: Analytically calculated running of scalar (to the left) and gaugino mass parameters (to 
the right), leading order only. The mass parameters are calculated as a function of M15 for the 
mSugra parameters m = 70 GeV and M 1/2 = 250 GeV. For M i5 ~ 2 x 10 16 GeV the mSugra 
values are recovered. Smaller Mi 5 lead to smaller soft masses in all cases. Note that the running is 
different for the different mass parameters with gaugino masses running faster than slepton mass 
parameters. 
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III. NUMERICAL RESULTS 



In this section we discuss our numerical results. All the plots shown below are based 
on the program packages SPheno 



591 and micrOMEGAs 



60, 



61 . We use SPheno V3 



including the RGEs for the 15 + 15 case [32|, [33| at the 2-loop level for gauge couplings and 



gaugino masses and at one-loop 
parameters, for a discussion see 



evel for the remaining MSSM parameters and the 15-plet 
3^| . For any given set of mSugra and 15-plet parameters 



SPheno calculates the supersymmetric particle spectrum at the electro-weak scale, which 



63j to calculate the relic density of the lightest 



is then interfaced with micrOMEGAs2.2 
neutralino, VL x oh? . 

For the standard model parameters we use the PDG 2008 values {jj, unless specified 
otherwise. As discussed below, especially important are the values (and errors) of the bottom 
and top quark masses, m& = 4.2 + 0.17 — 0.07 GeV and m t = 171.2 + 2.1 GeV. Note, the m t 
is understood to be the pole-mass and m b {m b ) is the MS mass. As the allowed range for 
^DAih 2 we always use the 3 o c.l. boundaries as given in 0|, i.e. Vt DM h 2 = [0.081,0.129]. 
Note, however that the use of 1 cr contours results in very similar plots, due to the small 
error bars. 

In the "seesaw sector" we have the parameters connected with the 15-plets, i.e. M 15 , Y 15 , 
Ai and A2. For the calculation of the dark matter abundance the most important parameter is 
M 15 . It has turned out that the effects of Y 15 , Ai and A 2 on the relic abundance of neutralinos 
are very minor. Note, however, that as discussed in the previous section, atmospheric 
neutrino oscillation data can not be explained in our setup, if the triplet mass is larger than 
approximately M 15 = Mt = 10 15 GeV. Also, the non-observation of lepton flavour violating 
(LFV) decays puts an upper bound on M 15 . The latter, however, is strongly dependent on 
tan/3 and depends also on m and M 1 / 2 . We will first show results using different values 
of M T as free parameter, without paying attention to neutrino masses and LFV. We will 
discuss how our results change for correctly fitted neutrino masses and angles towards the 
end of this section, where we also discuss and compare LFV excluded regions with DM 
allowed ones. 

We define our "standard choice" of mSugra parameters as tan/3 = 10, A = and fi > 
and use these values in all plots, unless specified otherwise. We then show our results in 
the plane of the remaining two free parameters, (m ,M 1 / 2 )- Fig. ((2j) shows in the top 
panel contours of equal dark matter density, fl x oh 2 . The lines are constant Q x oh 2 with 
Q x oh 2 = 0.1, 0.2, 0.5, 1,2. In the bottom panel we show the range of parameters allowed by 
the DM constraint at 3 a c.l. In both cases, to the left a pure mSugra calculation, whereas 
the plot to the right shows mSugra + 15-plet with Mt = 10 14 GeV. In each plot the yellow 
regions are eluded either by the lighter scalar tau being the LSP (to the bottom right) or 
by the LEP limit on the mass of the lighter chargino (to the left), m x + > 105 GeV. In 
addition, we show two lines of constant lightest Higgs boson mass, m h o = 110 GeV (dotted) 
and m h o = 114.4 GeV (dashed), as calculated by SPheno, see the discussion below. 

The plots show three of the different allowed regions discussed in the introduction. To 
the right the co-annihilation region, here the lightest neutralino and the lighter scalar tau 
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800 1000 



M 1/2 (GeV) M 1/2 (GeV) 

FIG. 2: Top: Contours of equal dark matter density (O x o/i 2 ) in the (mo, Mj^) plane for the 
"standard choice" tan (3 = 10, Aq = and /i > 0, for mSugra (left panel) and type-II seesaw 
with M T = 10 14 GeV (right panel). The lines are constant £\o/i 2 with ^ x oh 2 = 0.1,0.2,0.5,1,2. 
Bottom: Range of parameters allowed by the DM constraint at 3 a c.l. To the left: mSugra; to the 
right: Mt = 10 14 GeV. For a discussion see text. 



are nearly degenerate in mass. The line going nearly vertically upwards at constant Mi/ 2 
is the "focus point" line. The small region connecting the two lines are the remains of the 
bulk region, which has shrunk considerably due to the reduced error bars on Q DM h 2 after 
the most recent WMAP data The focus point line is excluded by the LEP constraint 
on the lighter chargino mass at low and moderate values of m . It becomes allowed only at 
values of mo larger than (very roughly) 1-1.5 TeV. However, note that the exact value of mo 
at which the focus point line becomes allowed is extremely sensitive to errors in m x + , both 
from the experimental bound and the error in the theoretical calculation. 

Comparing the results for the pure mSugra case to the mSugra+15-plet calculation, two 
differences are immediately visible in fig. ([2]). First, the focus point line is shifted towards 
larger values of M 1 / 2 - This is due to the fact that for the 15-plet at M 15 = 10 14 GeV the 
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neutralino is lighter than in the mSugra case at the same value of M1/2, compare to fig. 
(DQ). Maintaining the same relation between Mi and /1 as in the mSugra case requires a then 
a larger value of My 2 - Note that for the same reason the excluded region from the LEP 
bound on the chargino mass is larger than in the mSugra case. Second on finds that the 
co-annihilation line is shifted towards smaller values of m . The latter can be understood 
from fig. (J3j) . 
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FIG. 3: Allowed region for dark matter density (0.081 < S\o/i 2 < 0.129) in the (m ,M 1 / 2 ) plane 
for the "standard choice" tan/3 = 10, Aq = and /j > 0, for five values from Mr, Mt = 10 14 
GeV (red), to Mt = 10 16 GeV (cyan), to the left. To the right: Variation of the mass difference 
(top lines) and of £lh 2 (bottom lines), as a function of Mr for four different values of 
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mo: (cyan), 50 (magenta), 100 (blue) and 150 GeV (green) for one fixed value of Mi/ 2 = 800 
GeV. The yellow region corresponds to the experimentally allowed DM region. 



Fig. ([3]) shows the allowed region for the dark matter density in the (m , Mi/ 2 ) plane 
for our "standard choice" of other mSugra parameters for a number of different Mt (to the 
left). The plot shows how the co-annihilation line moves towards smaller values of mo for 
smaller values of Mt- The plot on the right in fig. ((3l) explains this behaviour. It shows 
the variation of the mass difference m fl — m xo (top lines) and of Qh 2 (bottom lines), as a 
function of M T for four different values of m : (cyan), 50 (magenta), 100 (blue) and 150 
GeV (green) for one fixed value of M1/2 = 800 GeV. The yellow region corresponds to the 
experimentally allowed DM region. Co-annihilation requires a small value of — rn Xo , 
typically smaller than a few GeV. With decreasing values of Mt the gaugino masses run 
down to smaller values faster than the slepton masses, thus effectively increasing m fl — m xo 
in these examples with respect to mSugra. To compensate for this effect at constant M 1 / 2 
smaller values of mo are required to get the m^ — m Xo in the required range. 

At this point a short discussion of the Higgs boson mass bound might be in order. LEP 
excluded a light Higgs boson with SM couplings with masses below m h < 114.4 GeV 3|. 
For reduced coupling of the Higgs boson to bb the bound is less severe, so this bound is 
not strictly valid in all of MSSM space. More important for us, however, is the theoretical 
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FIG. 4: Limits for mSugra with tan/3 = 10, and \i > for Aq = —300 GeV (left panel) and A$ = 
—500 GeV (right panel). The blue regions are allowed by the DM constraint, for the explanation 
of the bounds see fig. J2)) and text. 



uncertainty in the calculation of the lightest Higgs boson mass. SPheno calculates m^o 
at two-loop level using DR renormalization. Expected errors for this kind of calculation 
including a comparison of different public codes, have been discussed in 



64. 



64 . As discussed in 



65j even at the 2-loop level uncertainties in the calculation of rrufi can be of the order of 



66 , which calculates the 



3 — 5 GeV. In this context it is interesting to note that FeynHiggs 
higgs masses in a diagrammatic approach within the OS renormalization scheme tends to 
predict higgs masses which are systematically larger by 3 — 4 GeV, when compared with the 
DR calculation. We therefore showed in fig. ((2j) two lines of constant Higgs boson masses. 
The value of m h o = 114.4 GeV is taking the LEP bound at face value, while the lower 
value of myfi = 110 GeV estimates the parameter region which is excluded conservatively^ 
including the theoretical error. Since the lightest Higgs boson mass varies slowly with m 
and Mi/2, even a relatively tiny change in m h o of, say 1 GeV, shifts the extreme values of 
the excluded region by ~ 50 GeV in M\/2 (at small m ) and by ~ 150 GeV in m (at small 
M 1/2 ). 

Moreover, it is well known that the calculated Higgs boson masses are strongly dependent 
on the mixing in the stop sector and thus, indirectly, on the value of A . This is shown for 
the case of a pure mSugra calculation in fig. (jlj). Here we show two examples for the DM 
allowed region and the regions disfavoured by the Higgs boson mass bound at m h o = 114.4 
GeV and m^o = 110 GeV. Larger negative A Q leads to a less stringent constraint (for fi > 0). 
Note, that all of the bulk region becomes allowed at A = —500 GeV, once the theoretical 
uncertainty in the Higgs boson mass calculation is taken into account. We have checked for a 
few values of M T that for the case of mSugra+15 the resulting Higgs boson bounds are very 
similar. We thus do not repeat the corresponding plots here. Comparing the calculations 
shown in fig. (TJl) and the mSugra calculation in fig. (j2J) with each other, one finds that the 
DM allowed regions are actually affected very little by the choice of A . We have checked 
that this is also the case for mSugra + seesaw type-II. 

As mentioned above the uncertainty in the top mass is important for the calculation of 
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100 150 200 250 300 160 180 200 220 240 260 280 300 

M 1/2 (GeV) M 1/2 (GeV) 



FIG. 5: Logarithmically scaled zoom into the focus point region. In red the allowed region for 
0.081 < £lh 2 < 0.129 and in cyan the allowed region due the variation of m top = 171.2 ± 2.1 GeV. 
The left panel is for mSugra case and the right panel for = 10 15 GeV. The other parameters 
are taken at our "standard" values. 



the relic density. At low and moderate values of tan (3 the exact value of m t affects mainly 
the focus point region. As fig. ((2]) demonstrates near the focus point line the relic density 
changes very abruptly even for tiny changes of Mi/ 2 . This is because a comparatively small 
value of /i is required to get a sufficiently enhanced coupling of the neutralino to the Z Q 
boson. In mSugra the value of /i is determined from all other parameters by the condition of 
having correct electro-weak symmetry breaking (EWSB) and usually leads to M 1) M 2 <C /i. 
In the focus point region /i varies abruptly, points to the "left" of the focus point region 
are usually ruled out by the fact that EWSB can not be achieved. Since m t is the largest 
fermion mass, its exact value influences the value of /i required to achieve EWSB most. 
The change of /i with respect to a change of m t then can lead to a significant shift in the 
DM allowed region of parameter space. This is demonstrated in fig. ((5j), which shows a 
zoom into the focus point region for pure mSugra (to the left) and mSugra + 15 (to the 
right). The variation of the top mass shown corresponds to the current 1 a allowed range 
(s3|. The pure mSugra is especially sensitive to a change of m t . At large values of m the 
uncertainty in "fixing" M±/2 from the DM constraint can be larger than 100 GeV in the case 
of mSugra. Given this large uncertainty it would be impossible at present to distinguish the 
pure mSugra case from mSugra + seesaw, if the focus point region is the correct explanation 
of the observed DM. Note, however, that in the future the top mass will be measured more 
precisely. At the LHC one expects an uncertainty of 1-2 GeV at a linear collider m t 
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could be determined down to an uncertainyy of 100 MeV 

We now turn to a discussion of large tan (3. At large values of tan (3 the width of the CP- 
odd Higgs boson A becomes large, Ya ~ Ma tan 2 /9(m 2 +m 2 ), and a wide s-channel resonance 
occurs in the region m x o ~ Ma/2. The enhanced annihilation cross section reduces Q x oh 2 
to acceptable levels, the resulting region is known as the "higgs funnel" region. In fig. ((6l) 
we show the allowed range of parameters in the (m^M^) plane for one specific value of 
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500 1000 1500 2000 

M 1/2 (GeV) 



FIG. 6: Allowed region for dark matter density in the (mo,M 1 / 2 ) plane for Aq = 0, /x > and 
tan/? = 45, for (from top to bottom) M T = 5 X 10 13 GeV (red), M T = 10 14 (green) and M T = 10 15 
GeV (blue). 



tan/3 = 45 and three different values of M T . As demonstrated, the higgs funnel region 
is very sensitive to the choice of M T . It is fairly obvious that varying M T one can cover 
nearly all of the plane, even for fixed values of all other parameters. We have calculated 
the DM allowed region for various values of tan (3 and found that the funnel appears for all 
tan/3 > 40, approximately. 
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The strong dependence of the higgs funnel region on Mr unfortunately does not imply 
automatically that if large tan (3 is realized in nature one could get a very sensitive indirect 
"measurement" of the seesaw scale by determining (mo,Mi/ 2 ). The reason is that the higgs 
funnel is also very sensitive to the exact value of tan (3 and to the values (and errors) of 
the top and bottom quark mass. The latter is demonstrated in fig. (jZj), where we show the 
DM allowed range of parameters for a fixed choice of tan (3 and Mr varying to the left (to 
the right) m t (nib) within their current 1 a c.l. error band. The position of the funnel is 
especially sensitive to the exact value of m&. Comparing fig. (JJJ) with fig. (ED one can see 
that the uncertainty in m;, and m t currently severely limit any sensitivity one could get on 
M T . However, future determinations of and m t could improve the situation considerably. 
For future uncertainties in m t see the discussion above for the focus point region. For m b 
reference [69| estimates that m& could be fixed to 4.17 ± 0.05 GeV, which might even be 



improved to an accuracy of Am b ~ 16 MeV according to 



70 



1000 



> 
o 




100 



400 600 
M 1/2 (GeV) 



CD 

CD 



1000 



1000 




100 



400 600 
M 1/2 (GeV) 



1000 



FIG. 8: Allowed region for dark matter density in the (mcM^) plane for the "standard choice" 
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range for negligibly small neutrino Yukawa couplings (red) and Yt fitted to correctly explain solar 
and atmospheric neutrino data (blue lines). To the right: the DM allowed range of parameters for 
3 different values of A2, A2 = 0.5 (red), A2 = 0.75 (green) and A2 = 1 (blue). Note the logarithmic 
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All of the above figures have been calculated using fixed values for Ai and A2 and negligibly 
small Yukawa couplings Yt. This choice in general does not affect the calculation of the DM 
allowed regions much. However, a fully consistent calculation can not vary Mr, Yt and A 2 
independently, since this will lead to neutrino masses and angles outside the experimentally 
allowed ranges. Since Y T is diagonalized by the same matrix as the effective neutrino mass 
matrix, m„, see the previous section, the measured neutrino angles provide constraints on 
the relative size of the entries in Y T . The absolute size of Y T is then fixed for any fixed choice 
of A 2 and M T , once the neutrino spectrum is chosen to be hierarchical or quasi-degenerate. 
In the numerical calculation shown in fig. ((Hj) we have chosen neutrino masses to be of 
the normal hierarchical type and fitted the neutrino angles to exact tri-bimaximal (TBM) 
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values [15|, i.e. tan 2 #Atm = 1, tan 2 9 Q = 1/2 and sin 2 = 0. This has to be done in a 
simple iterative procedure, since the triplet parameters are defined at the high scale, whereas 
neutrino masses and angles are measured at low scale. For more details on the fit procedure 
see (i^ l . 

In fig. @ to the left we show two calculations of the DM allowed regions. The allowed 
range for negligibly small neutrino Yukawa couplings is shown by the filled (red) region, 
while the calculation with Yt fitted to correctly explain solar and atmospheric neutrino data 
is the one inside the (blue) lines. Note the logarithmic scale. As demonstrated, the exact 
values of Yt are of minor importance for the determination of the parameter region allowed 
by the DM constraint. Slightly larger differences between the fitted and unfitted calculations 
are found pushing M T to larger values (see, however, below). For smaller values of M T , the 
entries in Yt needed to correctly explain neutrino data are smaller and, thus, Yt affects the 
DM allowed region even less for M T < 10 14 GeV. 

In fig. (jHJ) to the right we compare three different calculations for A 2 , A 2 = 0.5 (red), 
A2 = 0.75 (green) and A2 = 1 (blue), for fixed choice of other parameters. This plot serves 
to show that also the exact choice of A 2 is of rather minor importance for the determination 
of the DM allowed region. Very similar results have been found for Ai, we therefore do not 
repeat plots ton$fflg WfYhfePfc 10 - A o=° < GeV > 
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FIG. 9: Allowed region for dark matter density in the (7710,-^1/2) plane for our "standard choice" 
of mSugra parameters and for two values of Mt- Mt = 5 x 10 13 (left panel) and for Mt = 10 14 
(right panel). Superimposed are the contour lines for the Br(fi — > cy). 



Finally, we will compare the constraints imposed on the parameter space of the model 
by Vtohih 2 with the constraints from the current data on non-observation of lepton flavour 
violating processes. Since LFV within the present model has been studied in some detail in 
331 ] . we will not repeat all of the discussion here. Instead, here we concentrate on fi —> 
exclusively, since the upper bound on Br(/j, — > cy) of £>r(/i — > cy) < 1.2 • 10 _n Q has been 
shown to provide currently the most important constraint. 

In fig. we show the DM allowed parameter regions for tan (3 = 10 and two values of 
M T , Mt = 5 ■ 10 13 GeV (to the left) and M T = 10 14 GeV (to the right), for a fixed choice 
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of all other parameters. Superimposed on this plot are lines of constant branching ratio for 
Br(n — > &y). The latter have been calculated requiring neutrino masses being hierarchical 
and fitted to solar and atmospheric neutrino mass squared differences and neutrino angles 
fitted to TBM values. Within the (m , M1/2) region shown, Br(fi — > cy) can vary by 
two orders of magnitude, depending on the exact combination of (m ,Mi/ 2 ), even for all 
other parameters fixed. The most important parameter determining Br(f.i — > ery), once 
neutrino data is fixed, however, is Mr, as can be seen comparing the figure to the left with 
the plot on the right. While for My = 10 14 GeV about "half" of the plane is ruled out 
by the non-observation of ji — > ej, for Mp = 5 ■ 10 13 GeV with the current upper limit 
nearly all of the plane becomes allowed. The strong dependence of \x — > &y on M T can 
be understood from the analytical formulas presented in (2^. In this paper it was shown 
that Br(fi — > ej) scales very roughly as Br{ji — >• ej) oc Mylog(Mr), if neutrino masses 
are to be explained correctly. For tan/3 = 10 one thus concludes that with present data 
values of M T larger than (few) 10 13 GeV - (few) 10 14 GeV are excluded by Br(/i — > ej), to 
be compared with < 10 15 GeV from the measured neutrino masses. Note, however, 

that (i) the constraint from neutrino masses is relatively independent of tan (3, m and M1/2, 
while fi — > e7 shows strong dependence on these parameters; and (ii) allowing the value of 
the reactor angle sin 2 9 R to vary up to its experimental upper limit, sin # R = 0.056 
leads to larger values of Br(fi — > ej) and thus to a tighter upper limit on M T . 
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In fig. ffTOl) we show the results for a calculation comparing dark matter and LFV in the 
case of large tan j3. Here the same constraints as in fig. ([9]) are shown, however for tan (3 = 45. 
Again we show the calculation for two values of Mr, since Mt is the most important free 
parameter. It is known that at large values of tan/3, LFV decays are enhanced due to 
an enhanced chargino diagram, which in the limit of large tan/3 scales approximately as 
tan 2 /3 2a] . Therefore, constraints on the parameter space from non-observation of LFV 
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decays are more severe in case of large tan/?, leading to tighter upper limits on Mt- This 
is clear if we compare fig. ^ and fig. (fTO]) . noticing the different scales. However, because 
of the higgs funnel region developing for large tan /3, the interesting part of the parameter 
space enlarges compensating for the larger values of the LFV decays. This can be seen in 
fig. (flOl l. where for M T = 5 x 10 13 GeV (left), most of the (m ,Mi/ 2 ) plane is allowed by 
the upper limit on Br{[i — > ej), while for Mt = 10 14 GeV (right), about "half" of the plane 
is ruled out by this limit. 

We have concentrated in this paper on discussing DM in mSugra with a seesaw type-II. 
Before closing this section, we would like to briefly comment on the case of seesaw type-I. In 
seesaw type-I one adds two or more singlet superfields to the superpotential of the MSSM. 
These singlets have Yukawa couplings to the standard model lepton doublet and a Majorana 
mass term, but no other couplings to any of the MSSM fields. The running of the mSugra 
soft parameters in this setup is therefore only changed by the neutrino Yukawa couplings. 
Just as in the seesaw type-II one can estimate from current neutrino data that the Yukawa 
couplings are order Y v ~ 0(1) for the right-handed Majorana mass order C(10 15 ) GeV. 3 
For any Mm smaller than this number, one therefore expects that the running of the soft 
parameters is essentially mSugra-like. (Apart from small off-diagonal terms in m|, which 
are exactly zero in mSugra.) This implies that also the DM regions should be very close 
to those found in the mSugra case. We have confirmed this expectation by calculating the 
DM allowed region for our standard choice of mSugra parameters and various values of the 
right-handed neutrino masses. Even for Y v at the upper limit allowed by perturbativity 
we did not find any significant departure from the mSugra case. With the hindsight of the 
results shown in fig. (JHj) for the seesaw type-II this is not surprising. 



One exceptional case for the seesaw type-I has been discussed, however, recently in [56 
The authors' of H observed that for Yukawa couplings dose to one and a large value of 
the common trilinear A , say A = 1100 GeV, the left sneutrinos can be the next-to-LSP 
(NLSP) for small-to-moderate values of mo, My 2 and tan/3. For a sneutrino NLSP nearly 
degenerate with the lightest neutralino a new co-annihilation regions then shows up at small 



values of mo. We have repeated this calculation with three right-handed neutrinos ([56J] use 
only one singlet superfield) and confirm the sneutrino co-annihilation region for \Y U \ ~ 0(1) 
and large A . However, in our calculation, if we insist on fitting the large atmospheric 
and solar angles, all of the region is excluded by upper limits on LFV decays, if we gut 
the matrix R of the Casas-Ibarra parameterization for the neutrino Yukawa couplings 



to the identity matrix. As has been shown in [27], |29|], we could, in principle, avoid these 
strong constraints from LFV by a careful adjustments of the unknown parameters in R. We 
did, however, not attempt to do a systematic study as to how R has to be chosen that the 
sneutrino co-annihilation becomes consistent with LFV decays. 



3 



Different from the seesaw type-II, where Y T depends linearly on M T , however, in seesaw type-I Yukawas 
scale like Y v ~ \[Mm. 
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IV. CONCLUSIONS 



In conclusion, we have calculated the neutralino relic density in a supersymmetric model 
with mSugra boundary conditions including a type-II seesaw mechanism to explain current 
neutrino data. We have discussed how the allowed ranges in mSugra parameter space change 
as a function of the seesaw scale. The stau co-annihilation region is shifted towards smaller 
mo for smaller values of the triplet mass Mt, while the bulk region and the focus point line 
are shifted towards larger values of Mi/ 2 for Mt sufficiently below the GUT scale. The higgs 
funnel, which appears at large values of tan f3 has turned out to be especially sensitive to the 
value of Mt- Determining iWi/ 2 from the mass of any gaugino and from a sparticle which 
is not important for the DM calculation, one could, therefore, get a constraint on Mt from 
the requirement that the observed flDAih 2 is correctly explained by the calculated Q x oh 2 . 

On the positive side, we can remark that current data on neutrino masses put an upper 
bound on M T of the order of 0(1O 15 ) GeV. Since this is at least one order of magnitude 
smaller than the GUT scale, the characteristic shifts in the DM regions are necessarily non- 
zero if our setup is the correct explanation of the observed neutrino oscillation data. Even 
more stringent upper limits on Mt follow, in principle, from the non-observation of LFV 
decays. A smaller Mt implies larger shifts of the DM region. However, the "exact" upper 
limit on M T from LFV decays depends strongly on tan/5, m and Mi/ 2 , and thus can be 
quantified only once at least some information on these parameters is available. 

On the down side, we need to add a word of caution. We have found that the DM 
calculation suffers from a number of uncertainties, even if we assume the soft masses to be 
perfectly known. The most important SM parameters turn out to be the bottom and the 
top quark mass. The focus point line depends extremely sensitively on the exact value of 
the top mass, the higgs funnel shows a strong sensitivity on both, m b and m t . 

Finally, it is clear that quite accurate sparticle mass measurements will be necessary, 
before any quantitative conclusions can be taken from the effects we have discussed. Unfor- 
tunately, such accurate mass measurements might be very difficult to come by for different 
reasons. In the focus point region all scalars will be heavy, leading to small production cross 
section at the LHC. In the co- annihilation line with a nearly degenerate stau and a neu- 
tralino, the stau decays produce very soft taus, which are hard for the LHC to measure. And 
the higgs funnel extends, depending on tan/3 and M T , to very large values of (m , Mi/ 2 ), at 
least partially outside the LHC reach. Nevertheless, DM provides in principle an interesting 
constraint on the (supersymmetric) seesaw explanation of neutrino masses, if seesaw type-II 
is realized in nature, a fact which to our knowledge has not been discussed before in the 
literature. 
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